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ABSTRACT 

The problem of synchronization and detection of random pulse-position-modulation 
( P P Ni ) sccj ucnccs js investigated under the assumption of perfect slot synchronization. 
Maximum-likelihood PPM symbol synchronization and receiver algorithms are derived 
that make decisions based both on soft as well as hard data: these algorithms are seen 
to be easily implementable. We derive bounds on the symbol error probability as well 
as the probability of false synchronization that indicate the existence of a rather severe 
performance floor, which can easily be the limiting factor in the overall system performance. 
The performance floor is inherent in the PPM format and random data and becomes more 
serious as the PPM alphabet size Q is increased. A way to eliminate the performance floor 
is suggested by inserting “special” PPM symbols in the random data stream. 
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I. INTRODUCTION 


Pulse-position-modulation is a modulation format known to be optimal in various 
ways for the direct-detection optical channel (see for example [1,2,3]). Under Q-ary PPM, 
information is contained in the position of a signal pulse in only one of Q subintervals, 
known as slots, dividing the symbol interval. 

When pulse-position-modulation is used in communication systems, the practice is to 
first achieve slot synchronization before attempting higher order synchronization and sym- 
bol decoding. Slot synchronization is usually obtained by using a tracking-loop as recently 
studied, for example, by Chen and Gardner [4] and Ling and Gagliardi [5j. Although op- 
timal (in the maximum-likelihood as well as the mean-square-error sense) PPM slot and 
symbol synchronizers have also been recently derived [6.7], they have the disadvantage of 
being more complicated to implement compared to the tracking-loop synchronizers that 
are well understood and easy to implement. The difficulty in implementing the optimal 
synchronizers st ems from the need to record the exact arrival time of each del ected photon, 
a task that may be difficult to achieve at high data rates and large signal intensities. A 
further reason that makes tracking-loops more desirable for slot synchronization is that 
they have been shown to result in receivers with symbol error probability performance 
within a fraction of a dB from the perfect slot synchronization case at reasonable signal 
levels [4]. We point out here that the authors in [4] are investigating the effects of slot 
synchronization errors only and thus assume that once slot synchronization is achieved, 
symbol synchronization is automatically obtained. This is equivalent to assuming that the 
only ambiguity in symbol synchronization is the ambiguity in slot synchronization, which 
in practice is not a valid assumption since slot synchronization does not imply symbol 
synchronization. 

In this paper we investigate the problem of PPM symbol synchronization and decoding 
under t he assumption of perfect slot synchronization. For Q-arv PPM, the existence 
of slot synchronization still leaves unresolved a Q-ary ambiguity as to the position of 
the PPM symbols. However, in contrast to the slot synchronization case where system 
error probability degrades gracefully with the slot synchronization error, the effect of non- 



perfect symbol synchronization is catastrophic. One can easily see that a sequence of N 
random PPM symbols decoded under non-perfect symbol synchronization will result in all 
N symbols being decoded erroneously, irrespective of the size of the error. This observation 
leads to the conclusion that the real bottleneck in system performance is due to the symbol 
synchronization subsystem which we investigate in the sequel. 

In section II we investigate and characterize the synchronization properties of random 
PPM sequences. In section III, we derive ML symbol synchronization algorithms for the 
PPM, optical Poisson channel, both from soft as well as hard data. Here we also derive a 
bound on the synchronization probability, valid at high signal-to-noise levels. Section IV 
contains the derivation of optimal receivers that make sequence decisions in the absence 
of symbol synchronization. Also included in this section is a bound on the minimum 
achievable symbol error probability. Finally, we conclude with section V. 

II. SYNCHRONIZ ABILITY OF PPM SEQUENCES 

In this section we investigate the synchronization properties of random PPM se- 
quences. Broadly speaking, these are properties associated with our ability to identify 
uniquely (or not) the location of PPM symbols within a sequence of such symbols when 
only slot boundaries are known. Before we proceed further, we introduce some definitions 
and terminology to establish a common ground and facilitate later analysis. 

Definition 1: A binary sequence of Q digits is said to satisfy the PPM constraint and, 
thus, be a Q-ary PPM symbol only if exactly one out of the Q digits is a “one”; Q will be 
referred to as the PPM alphabet size. 

It is easily seen from the above definition that for a given sequence length N there are 

valid PPM sequences and that they can be thought of as a subset of the 2 N ® binary 
sequences of length NQ. 

Definition 2: A sequence of binary digits is said to be a valid Q-ary PPM sequence if. 
starting from the first bit. cv < ry consecutive subsequence in it of length Q is a Q-arv PPM 
symbol. 
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Definition 3: Consider a Q- ary PPM sequence (N + 1) symbols long ((N + 1)Q slots) 
and a sliding window of NQ slots applied to the PPM sequence. For a given slot shift j, 
j = 0, 1,2, — 1), we will refer to the binary sequence within the window as a PPM 

binary subsequence at (slot shift) j. It is obvious from the definition that a PPM binary 
subsequence is not necessarily a PPM sequence, i.e., not all consecutive binary sequences of 
length Q within the window are PPM symbols. Also obvious is the fact that if j — 0, then 
the PPM binary subsequence is a PPM sequence that coincides with the first N symbols 
of the original sequence. Similarly, for j = Q the PPM binary subsequence is a PPM 
sequence that coincides with the last N symbols in the original sequence. Notice that slot 
shift j = Q does not need to be investigated in searching for symbol locations since it is 
equivalent to j = 0. 

Given the above definitions, we are now ready to derive some results. As a first step in 
characterizing t he synchronization properties of random PPM sequences, we are interested 
in the following problem. 

Consider a Q-ary PPM sequence of length (A* 4- 1) symbols, chosen at random from 
the set of + possible PPM sequences, and a given slot shift j . We are interested in the 
probability that exactly K out of the N symbols within the PPM binary subsequence at 
slot shift j are Q-ary PPM symbols. The usefulness of this quantity should be apparent in 
the context of symbol synchronization when only slot synchronization is present. Leaving 
the details of the derivation for Appendix A, we can show that this probability, defined by 
P{K\ j,Q,N), is given by 

P(K-. j. Q, JV)= (£) {(j/Q) IK+ "(l - j/Q) tN ~ K> + (l - j/Q) <K+,, (i/0) <iv “ K| } • 

(ll 

It is obvious from the above expression that P(K\ j . Q, N ) = P(K; Q — j, Q, N), which 
implies that P(K: j , Q , N) is symmetric with respect to forward and backward slot 
shifts from the correct slot for symbol synchronization. Another observation that can be 
readily made is that the above probability is a convex function of j, achieving a minimum 
at j = Qi 2 when Q is even, as is usually the case. This implies that slots closer to the 
correct symbol synchronization slot have a higher probabilit y of being erroneously chosen 
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for symbol synchronization than slots further away. 

Of special interest is the case when K — N, P(N ; j, Q, N), i.e., the probability that 
the PPM binary subsequence at slot shift j is a PPM sequence. Easily obtained from (1), 
it is 

P(N-,j,Q , N) = (j/Q) N+ ' +(l-j/Q) N+ '. (2) 

Equation (2) implies that there is a nonzero probability of identifying symbol locations 
erroneously, even in the absence of noise. A bound on the probability of erroneous symbol 
synchronization is derived in the next section. 

In passing, we note that P{ K ; j. Q, N) can be expressed as the ratio of the number 
of PPM sequences that, for a given j, result in exactly K PPM symbols within the PPM 
binary subsequence at j , to the total number of distinct PPM sequences of length (N -4-1). 
Denoting the number of sequences resulting in K matches bv D(K; ?, 0, N), we then 
have 

D( A ; j. Q. N) = Q iN+ ' ] P(K : j. Q . N). (3) 

We now turn our attention to the implications of (1) on symbol synchronization in 
the limit as the sequence length N tends to infinity. It is readily seen that for any fixed 
K,P(K ; j , Q. N) — > 0 as N — > oo, which implies that the fraction of sequences that 
result in exactly K matches diminishes with N. However, the more important quantity is 
the probability that the number of matches K will exceed a given number r. Specifically, 
of interest is the smallest value of r such that the probabilit y of more than r matches goes 
to zero as .V —+ oo. for all slot shifts j = 1, 2. ■■■»((? — 1). Obviously, this minimum value 
of r is a function of the PPM modulation format and the alphabet size Q only; moreover, 
the smaller it is, the better the symbol synchronizer will be able to perform in the presence 
of noise for large N . 

The following proposotion establishes the range of values of r such that the probabilit y 
of more than r symbol matches tends to zero with N. 

Proposition 1: Let Pr[A" > r; j,Q, JV] be the probability that the number of symbol 
matches K within a PPM binary subsequence at slot shift j exceeds r for a given Q and 



A T . Then, if (r/N) > (1 - l/Q) 


lim Pr[K > r; j, Q, AT] = 0 (4) 

iX =? oc 

for all slot shifts j. Conversely, for (r/N) < (1 - l/Q), there is some slot shift j such that 

lim Pr\K > r; j, Q , N] = 1. (5) 

N=>oo 

Proof: We first prove (4). Using the Chernoff bound we have 

PrlK > r; j , Q , N\ < exp (—sr)<p{s; j, Q, N) s > 0. (6) 


where 


Q(s ; 




7\ r \ 
; 


iV 



TJJ — 0 


= 07Q)[0'/<?)‘'+d ~J/Q)f 


(i -;/Q)!d ->/(?)£•* -;/Q]' v (") 


is the characteristic function of the random variable A*, easily obtained from (1). In order 
to guarantee that (4) holds for all j. it is enough to make sure that (4) holds for the j 
that makes the right-hand side of (6) largest. From the easily established convexity of 
o(s; j , Q. N) and its symmetry around j = Qi 2, we conclude that the maximum occurs 
at one of the boundaries, say j — 1. Letting p = l/Q we have 

o(s; j , Q. JV) < </>($; 1, Q. N ) = p[p( a + (1 - p)] N + (1 —/>)[(! - p)e 3 4- p}* ■ (8) 


Observing now that the second term in (8) is the larger term we can further bound the 
right-hand side of (8) by twice this term; combining with the bound in (6), we get 


Prf A > r: j. Q. N] < 2(1 — p) exp( —sr) [(1 — p)e s — p]^ . s > 0 


(9) 


Further tightening the bound with respect to s, we finally obtain for (r/N) > (1 - l/Q) 

0 < Pr'K > r; j, Q. N] <2(1 -p)exp { - N[E(r/N\\( 1 - p)) - h b (r/N))} . (10) 
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In (10) i?[g||fi] and h b (x) are the cross entropy (inaccuracy) and binary entropy functions 
defined respectively by 


E[q\\h) = -q ln(h) - (1 - q)ln( 1 - h) (11) 

h b (x) = -x ln(x) - (1 - x) In (1 — x) (12) 

for some probabilities q , h and x. It is well known that the inaccuracy is always greater 
than or equal to the entropy (see for example [8]), which implies that the exponent in (10) 
is non-positive. Thus, in the region (r/N) > (1 — 1 /Q) where the bound is valid, taking 
the limit as N — > oo yields (4). 

We now turn our attention to the converse in (5). Since we only need to show con- 
vergence to unity for some j, we choose j = (# — 1) (an educated choice). We start by 
deriving the following C’hernofF bound 

PrlK < r; j, Q. N\ < exp( -sr)o{ _s; j , Q, N ), s < 0. (13) 

Following arguments paralleling those above we obtain for j = (Q — 1) and (r/N) < 

(i - i /Q) 

Pr\K < r; j, Q, .V] < 2(1 - l/Q)exp{-N[E(r/N\\( 1 - 1 /<?)) - h b (r/N)}} (14) 

where we have used the fact that Pr[K < r; j , Q , N] < Pr[K < r; j, Q, N}. The 
condition (r/N) < (1 - l/(?) gives a range of values of r for which the bound is valid 
and is derived from the condition .s < 0. Taking limits as N — > oo in (14) we obtain 
Pr[K < r; j , Q . .V] — > O'when (r/N) < (1 - 1 /Q). which in turn implies (5). This 
completes the proof of the proposition. 

One of the implications of the theorem is that for large (theoretically infinite) A’, the 
fraction of symbols within any PPM biliary subsequence that are valid PPM symbols can 
be quaranteed to be less than (1 - 1 /Q). In other words, the number of PPM symbol 
matches K for any slot shift j = 1, 2, ••*,((? — 1). satisfies, for large N 

K<(1-1/Q)N (15) 
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which implies that the number of mismatched symbols, call it d , satisfies 

d = (N - K) > N/Q. (16) 

From (16) we can see that for large A T , the minimum number of symbol mismatches for 
any slot shift j , call it d m ; n , is 

4n„ = N/Q. (17) 

The quantity d m in can easily be paralleled to the minimum distance of a block code, only 
now d m in is the minimum average distance between random PPM sequences for any slot 
shift j = 1. 2, •••,(<? — 1) away from the correct synchronization slot. Clearly, the larger 
d m in is. the better the symbol synchronizer will be able to identify symbol locations in the 
presence of noise. Another obvious observation from (17) is that d m ; n is monotonically 
decreasing with the PPM alphabet size Q, which in turn implies that synchronization 
performance deteriorates as Q is increased. This observation was made also in [4] and [7] 
for the slot and PPM symbol synchronization problems respectively. Later in this paper we 
relate the minimum distance d m j n for random PPM sequences to the smallest achievable 
synchronization error probability. 

In concluding this section, we note that the minimum average distance d mm derived 
above is for random PPM sequences, which indicates the existence of sequences with larger 
dmin- The problem of designing such sequences is briefly investigated in the next section. 

We next turn our attention to the problem of deriving optimal PPM symbol synchro- 
nization when slot synchronization is present. 

III. ML SYNCHRONIZATION 

Our application area here is the direct-detection optical Poisson channel. For this 
channel, it can be easily established that when only slot synchronization is present, the 
sufficient statistic is the number of photons (events) observed in each slot interval of T' 
seconds. For the Poisson direct-detection channel, the number of counts observed in dis- 
tinct slot intervals are independent Poisson random variables with intensity (A* + A„) if a 
signal pulse is present in the corresponding slot and A n otherwise: A s and A n are known 
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as the signal and noise intensities respectively. We derive synchronization algorithms from 
two kinds of observations, defined in the sequel. 

a) Soft-data observations: 

In deriving the ML synchronization algorithm here, we assume that the receiver has 
available an observation vector R = (A'i, A*2. •••, Kq, A 'nq) with elements being 

the number of counts in NQ observed, consecutive slots. We will refer to the data in R as 
soft-data. Within this observation window there are (N + 1) symbols, two of them at the 
boundaries partially contained. The problem is to estimate the location of symbols within 
the observation interval. 

Denoting by m the variable associated with the start of a symbol within the first Q 
slots in the observation interval, a ML synchronizer implements the following 

max V' Pr( R/ m. d) (18) 

0<m<(Q-l) - — ' 

d 

where d = (rfj. d 2 . dQ. ,di\ - q) is the random modulation vector of intensities 

within the observation interval. In writing (18), we assume that all modulation vectors 
are equiprobable. Following the general approach described in [6, 7] and making the 
assumption that the first, partially contained symbol within the observation interval is a 
continuation of the last partially contained symbol, we obtain as the likelihood statistic 

L s (m) = J^ In jv ex P!^(^)-^(. — DQ+j+mlj • (19) 

In (19). A* = (1 — As /'A n ) and K, is the number of counts observed in the i - ih slot in the 
observation interval. The indexes in (19) are interpreted modulo QN to account for our 
approximation above. We note here that this approximation was not necessary in deriving 
the ML synchronizer but was made to reduce complexity at no practical performance loss 
for values of of N greater than four [6. 7j. 

A further approximation to (19) can be obtained by using only the largest term in the 
sum over j to yield 

.v 

L (m) = A( l _i)Q_ rJ ;._|_ m . (20) 
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We will refer to (20) in the sequel as the max-rule. In (20), ji is the value of j, 1 <j<Q, 
that, for a given i and m maximizes K^^DQ + j +m . 

Computer simulation results comparing the probability of correct synchronization, 
P cs , for the synchronizers in (19) and (20) are reported in Figures 1,2 and 3 for different 
values of Q and N. These results indicate that the much simpler synchronizer in (20) 
performs practically as well as the significantly more complicated one in (19). Also ev- 
ident from the graphs is that for the same signal energy per slot \ S T' , synchronization 
performance degrades with Q and improves with N , as predicted in the previous section. 
Finally, it is clear from the figures that a performance floor exists which is practically in 
effect for signal levels of about A ,T' = 10.0. This error-floor is investigated later in this 
section. 

We now turn our attention to deriving synchronization algorithms when the receiver 
observations consist of hard-data, obtained from the soft-data vector R by making hard 
decisions in each slot interval. 


b) Hard-data observations 


Here we assume that the receiver has available the binary vector of observations X = 
(a*!, x 2 , xq, •••, xnq) obtained by making hard-decisions on the vector R. By 

this we mean that the receiver looks at A',, 1 = 1, 2, •••, NQ. and decides that x x = 1 
or x , = 0 according to whether a signal pulse is detected (a “one”) or not (a “zero”) 
respectively; i.e.. 

Pr f “one” /i?, = A',1 * i = 1 

— — — > 1 ( 21 ) 

Pr\*zeTO*/Ri = Ki] K 1 1 

In (21). R, is the (Poisson) random variable associateed with the number of counts in the 
i-th slot and A\ are the observed counts. The above test can easily be shown to reduce to 

^ i < 1 ( 22 ) 

* t =u 

where the optimal threshold -y is given by 


A, r -r In (Q-l) 
ln( 1 — Xs/\n ) 


10 


It is obvious that further processing on the vector X to derive synchronization and perform 
decoding is much easier to implement compared to the soft-data case. Another advantage of 
algorithms derived from hard-data is that their structure does not depend on the complete 
statistics of the channel. The above reasons were partly why hard-data was employed in 

i»;- 

Before we proceed with the derivation of the ML synchronization algorithm, we intro- 
duce the following quantities: 

Pr\Xi = 0/di = 0] = Pr[Ki < T/di = 0] = P 0 „ (24) 

PrLY, = 1/di = 1] = Pr\K, > 7 /d, = 1] = P n . 

Then 

Pr\Xi-\l di = 0] = (1 — Poo) (25) 

Pr X, = 0/dj = 1 = (1 - Pn )• 


The above probabilities can be easily precomputed given T' , Q and the signal and 
noise intensities A, and A„ respectively. Under the assumption of equiprobable PPM 
symbols, the ML synchronization rule maximizes over 0 < m < (Q — 1) the following 
statistic 


L 1 ( 77 ? ) = Z *•[*/ m, d] 

d 

N-l Q Q 

— || ^ P r [^tQ+ ] + rn ~ ZiQ + j + rn / W* d t Q~ j + m — 0 ] 

t = u A = 1 ; = i 

X Pr[A t Q + k + m — •TfQ + fr-f-m d iQ+ k + m — lj (26) 


N-l Q 

— 1 1 PT[AtQ->-j+rn — 'EiQ-rj-hm/m* + = 0 : 

i=0 j = 1 


Pr X,Q + k + m — 'I'iQ^k + m / djQ + k-r-rv — lj 

Pr X { Q^k± m — iQ+ JL’ttti / rn . * = 0] 


A = 1 


The second equality above is a result of the independence between choice of PPM sym- 
bols. which breaks the expectation over sequences to expectations over individual symbols. 
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The third equality should be self evident. Observing now that the product over all j ' s and 
is of Pr\XiQ+j+ m = x t Q+j +m /m, diQ + j +rn = 0] is not a function of m, we can equiva- 
lently maximize, 


N - 1 


L 2 (m) = £ 


1-0 


£ '! 


iQ+fc + m 


U= 1 


where 




Pc[-\»Q-|-Jfc-t-m — %iQ+k + m /m , + m — 1] 


(27) 


(28) 


iQ+k+m Pr[XiQ+k+m = Zj <?+ fc + m /”*> d ; Q + fc +m = 0] 

In deriving (27), we have dropped terms not dependent on m and taken the logarithm 
of the resulting expression. Further noting that by multiplying f'g +fc+m by any constant 
does not affect the maximization, we finally obtain 


A' — 1 


Lh{m) — y; In 


1=0 


,k—\ 


(29) 


where 


Pn 




U Q+fc + m - D l iQ + ; + ni - 

1 1 ~ -COO ) 

It is clear that l,Q+k- i-m takes one of two values according to 


(23) 


h = 


1, if a-, = 1 

C, if x , = 0, 


where 


c = 


(1-P„)(1 - Poo) 


(31) 


(32) 


Pn Poo 

Since for a practical system Pn > (1 — Pn) and Poo > (1 — Pooh it is 0 < C < 1. For a 
good system (Pn >> (1 — Pn ) and Poo >> (1 — Poo)) C will be very close to zero and 
exactly zero when either Poo or Pn is unity. In this case, the optimal synchronizer from 
hard data reduces to 


Lh(m ) 


A’-l 

£ 

/ = 0 


In 


Q 

v 

U = i 


X iQ + k + m 


(33) 


A situation when 0 = 1 arises, for example, when = 0 resulting in an erasure 
channel. It is then seen that the optimal synchronization rules for a perfect channel where 
both ^zeros’’ and "ones" are decoded correctly with probability one is the same as that 
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for the erasure channel where only “zeros” are detected correctly all the time. This of 
course does not imply that the performance of the synchronizers for the two channels are 
identical. From the information theory point of view, for both the perfect and the erasure 
channels, it is the presence of a pulse that carries all synchronization information, utilized 
fully by (33). 

In general, for channels with small noise intensities A n and relatively large signal 
intensities A a , (33) can be used as an approximation to (29) to further reduce complexity. 

Having derived optimal synchronization under the assumptions of both soft as well 
as hard-data, we now turn our attention to the investigation of the ultimate performance 
achievable by such synchronizers. 

c ) A bound on synchronization probability 

Our interest here is deriving an upper bound to the probability of correct synchro- 
nization. P, . It is clear that the synchronization probability is bounded from above by 
the probability that one or more PPM binary subsequences are PPM sequences, in which 
case a random selection must be made. Denoting this random sequence limited probability 
by Prsi - we have 

Pc»<p r ,i = Q- {N+1) J2 r^T’ (34) 

fc=o fc + 1 

where P*. is the number of PPM sequences resulting in exactly k slot shifts for which PPM 
binary subsequences are valid PPM sequences. The probability of k such matches, P m (fc), 

is 

P m (k) = D k (35) 

Leaving the details for Appendix B. we can derive the following expressions 

D 0 = <? (A ' +1) - 2 (Q - l) (A '~ n - (Q - 2) ( A + 1) (36) 

P m (0) = 1-2(1 - 1 Q) , V ' H - (1 - 2/Q) ( ' v+1) (37) 


Dl ~ { 2[((? - lp' +1 > - (Q - 2) 


.v+n 


-11 


Q = 2 
Q> 2 


(38) 


13 



The derivation of the above equations was obtained by repeated use of the following 
observation. 

Proposition 2: For a given PPM sequence, in order for the PPM binary subsequence at 
j to be a valid PPM sequence it is necessary and sufficient that either all pulses (pulsed 
slots) are before the j-th slot or all after the _/-th slot. 

The proof of the above proposition is easily seen by construction. 

YVe are now ready to derive a bound on the random sequence limited probability P r »i, 
and by extension to the probability of correct synchronization P ca . 

Proposition 3: The following bound on the random sequence limited probability P ra i 

holds for all Q and N 

Prsi<l-d- 1/Q) (N+1) + ^(1 - 2/Q) (N + 1) (42) 

with equality for Q = 2 and Q = 3. 

Proof: The following inequality is obvious 

Q-2 

Prsl < P m( 0) + \ P m { 1 ) ~ ^ PM ~ 1) + l ~ Y PM) (43) 

where equality is when Q equals two or three. Using the fact that 

Q-2 

V P w (A-) = l-|Pn,(0)-rP IM (l)-P nl «?-l)] (44) 

k = •> 

and t he expressions for the various probabilities in equations (37). (39) and (41 ) we obtain 
the required bound. 
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We note that the bound in (42) will be tight for all values of Q when N is sufficiently 
large to make the probability of two or more slot shifts for which PPM binary subsequences 
are PPM sequences negligibly small. This fact is verified in Table 1, where simulation 
results are compared to our bound. 


A lower bound on the probability of erroneous synchronization P e , = (1 — P cs ) is 

Pea > (1 - 1 /Q) iN + 1) ~ ^(1 - 2/C?) (N+1) (45) 


which is achieved at high signal- to-noise levels and sufficiently large N. It is clear from (45) 
that the error-floor rises with increasing Q and decreasing N as predicted by the decrease 
in d min in equation (17). Let us now show that the distance d m j n enters the bound on 
erroneous synchronization probability in a natural way as N — ♦ x. For a given d m i„, we 
have Q = N/d min . Then 


lim P ( . , > lim 

_\ — oc .V — oc 


rfm 1 /, 2 d min \ S *'\ 

-—) 2 \ 1 ~X ) / 


— t 


-d n 


\ £ 2 <i m in 

3 


(46) 


Since the second term above is at least an order of magnitude smaller than the first for 
reasonably large d min ( say d min > 3), we obtain, for very large (theoretically infinite N) 
and reasonably large d m i n 

Pe, > e~ dtl " n . (47) 


Equation (47) indicates an exponential decrease in the error-floor with increasing d m j n . 

As an example of how the bound in (45) can be used in a system design we derive 
next the minimum number of slots that need to be processed in order to quarantee a given 
performance. The number of slots observed, W = Q{ AT + 1 ), is a measure of the receiver 
complexity and is an important parameter in practical implementation. We show that the 
following proposition is true. 

Proposition 4 : The smallest number of slots ir min that need to be processed to guarantee 
an error-floor not more than some P e f is 


WVnin = “2 log 2 (P e /). 


(48) 
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and is achieved with binary PPM. 

Proof: We have W = Q(iV + l)>2(iV + l) with equality iff Q = 2. Further, we have 

(1 - 1 /Q) iN+1) - ^(1 - 2/<?) (N+1) = P ef 

which implies (A" + 1) > - log 2 (P e /) with equality iff Q = 2. Combining the two inequal- 
ities above, we get 

> —2 log 2 (P e/ ) 

with equahty iff Q = 2. This proves the proposition. 

We end this subsection bv noting that the bound in (42) is valid also for the re- 
sults reported in [12] wdiere PPM symbol synchronization is derived in the absence of slot 
synchronization. Next, we investigate ways to remove the error floor predicted above. 

d) Sequence Design 

We start with a definition. 

Definition 4: A PPM sequence will be said to be svnchronizable if no PPM binary sub- 
sequence in it is a valid PPM sequence. Equivalently, a PPM sequence is synchronizable 
if each PPM binary subsequence of it contains at least one symbol that is not a PPM 
symbol. 

It is clear from the above definition that symbol locations within synchronizable PPM 
sequences can be uniquely identified in the absence of noise. As we derived earlier, for a 
given Q and A' there are D 0 synchronizable PPM sequences as given by equation (36). 

We are now ready to prove the following propositions. 

Proposition 5: When a synchronizable PPM sequence is inserted in another PPM se- 
quence, the resulting longer sequence is synchronizable. 

P roof: By the definition of a PPM sequence all symbols within it must be PPM symbols. 
However, for any slot shift j = 1. 2. ••• , (Q - 1). there is at least one symbol in the 
inserted synchronizable sequence that is not a PPM symbol (by definition). Since the 
inserted sequence is a part of the longer sequence, this implies that at least one symbol in 
the longer sequence is not a PPM symbol for all j = 1. 2, ••• , (Q — 1). This implies by 
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definition 4 that the longer sequence is synchronizable, which completes the proof. 
Proposition 6: A (?-ary sequence is synchronizable if and only if it contains both symbols 
1 and Q ; symbol 1 has the pulse in the first slot and symbol Q has the pulse in the Q-th 
slot . 

Proof: a) We first prove the forward statement. If both symbols 1 and Q are in the 
sequence, then there is no slot shift j = 1, 2, • • ■ , (£? — 1) for which all pulses are before 
or after j. Thus, by proposition 2 no PPM binary subsequence is a PPM sequence and, 
thus, by definition 4 the sequence is synchronizable. 

b) If not both symbols 1 and Q are in the sequence, we distinguish three possibilities: 
only symbol 1 or only symbol Q is present or neither is present. If symbol 1 only is present, 
then the statement "all pulses are before slot ((? — 1)’“ is true which by proposition 2 it 
implies the sequence is not synchronizable. Similarly, if symbol Q only is present, then 
the statement “all pulses are after slot 1" is true which implies again that the sequence 
is not synchronizable. Finally, if neither symbol is present, both of the above statements 
in quotation marks are true, which again implies the sequence is not synchronizable. This 
completes the proof. 

To facilitate reference in the sequel, we will refer to the pair of symbols 1 and Q as 
the synchronizable pair . 

It is clear from the above proposition that Do in equation (36) gives the number of 
sequences from the set of possible sequences that contain at least one synchro- 

nizable pair. Obvious from the proof of the proposition, also, is the fact that it doesn't 
matter where the two symbpls are located within a sequence to make it synchronizable. 
Finally, we observe that the more synchronizable pairs a sequence contains, the better its 
synchronization properties in the presence of noise. 

One way of removing the error-floor predicted by (42) is to periodically insert in 
the random data stream a synchronizable pair. If. for example, a synchronizable pair is 
inserted every L PPM symbols, all PPM sequences of length greater or equal to L will be 
synchronizable according to proposition 5. The efficiency E of such a scheme as measured 


by the number of information symbols per transmitted symbol will be 


E=j±^=l-2/L, ( 49 ) 

which approaches unity as L is increased. It is clear, however, that the smaller L is, 
the better chances for correct synchronization in the presence of noise will be. In cases 
where a special synchronization pattern is inserted in the data stream to facilitate frame 
synchronization (see for example [10, 11, 12]), the extra reduction in efficiency to improve 
symbol synchronization can be avoided by choosing a pattern with as many synchronizable 
pairs as possible. Such a pattern will aid in both symbol as well as frame synchronization 
when they are obtained separately. 

In the next section, we investigate the problem of deciding what the symbols within 
an observation interval are as opposed to where they are which we analyzed above. We do 
this under the assumption of slot synchronization only. 

IV. DETECTABILITY OF PPM SEQUENCES 


a) ML Receivers: 

In this section we turn our attention to the problem of PPM sequence estimation 
when only slot locations are exactly known. Our observations consist of the vector R = 
( A'i . A '2 , ■ • • , Kq . • • • , A'/vq ) of slot counts. 

Under the assumption of equiprobable sequences, an optimal receiver is a ML receiver 
that chooses as its sequence estimate the sequence d that maximizes 

Q - 1 

p(R/d) = p(R/d, m )Pr(m ). (50) 

m — 0 

In (50). Pr(m) is the a priori probability that the first PPM symbol in the interval [0.QT' 
starts at time ml'. Assuming that no a priori knowledge exists, then Pr(m) = 1/Q for all 
m. Using the fact that the elements of R are conditionally independent Poisson random 
variables and dropping terms not dependent on the modulation sequence, we obtain as our 


optimal receiver 

Q-i r N 

max £(d)=Vexp ln( X) V A'j. (d)+m . (51) 

deA L ^ J 

In deriving (51) we have assumed that the first symbol in the observation interval is a 
continuation of the last, partially contained symbol. As for the synchronization case, this 
approximation reduces substantially the complexity of the receiver at no practical per- 
formance loss for N greater than about four [ 11 ]. To accommodate this approximation, 
indexes in (51) are interpreted modulo QN. The set A is the set of all possible distinct 
Q-ary PPM sequences of length N and J,( d), i = 1 , 2 , • • • ,N, is a set of indexes indicat- 
ing the slots within sequence d that contain the pulses. For example, the Q = 4 PPM 
sequence consisting of the N = 5 symbols {1.3. 3. 4. 2} can be described equivalently by 
{J,(d)} ( A j ={1.7.11.16,22}. Finally A’ = (1 + \ a /\ n ) is as defined earlier. 

Although each computation of (51) is relatively easy to perform, the complexity of 
the receiver is still overwhelming since statistic need to be evaluated before a decision 
is made. However, following arguments similar to those in 111:, we can show that only Q 
of the Q N sequences are most likely to have been sent. Denoting the set of Q candidate 
sequences by A”, the receiver in (51) becomes (with no loss in performance) 

max f'(d) . (52) 

dPA* 

The set of Q candidate sequences A* is obtain in the following way: 

Step 1: Group the observation vector into Q consecutive slots per symbol and decode it into 
a sequence of N PPM symbols. Store the decoded sequence as a possible candidate. 
Note that decoding the observation vector into PPM symbols is done by choosing the 
largest number of counts in each group of Q slots, which is the optimal strategy when 
synchronization is present. 

Step 2: Cyclically shift the observation vector by one slot to the left and go back to step 1. 
Repeat until the observation vector is cyclically shifted by {Q — 1) slots. 

The result of the above procedure is Q candidate sequences which are then used 
to evaluate f'(d) and make a final decision by choosing the largest. The idea behind the 
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tremendous reduction in the number of evaluations of £(d) is simple: we know that symbols 
start at one of the times mT' ,m = 0, 1 • • • ,(Q — 1) corresponding to the start of the first 
Q slots in the observation interval. Thus, by decoding the observation vector for each 
possible 77? = 0, 1, • • • , (Q — 1 ), we are assured that one of the Q decoded sequences is 
the one we would have obtained have we had perfect synchronization. This implies that 
limiting our search in the set A“ does not increase our error probability, since the best we 
can hope to achieve is the performance of a perfectly synchronized receiver. 

To illustrate the above approach, we take a simple example. 

Example 1: Let Q = 3, N = 4 and the observation vector R = (0. 4, 2, 1, 5, 3, 2, 3, 1, 6, 3, 5, ) 
Decoding this vector into a PPM sequence, we obtain do = (2,5,8,10). Cyclically shifting 
R by one slot to the left and decoding we obtain dj = (1,4,9,11). Finally shifting by 
another slot we get d 2 = (3,4,8,10). Evaluating (52) using ln(X -) = 1.3 we get £(do) = 
1.462 * 10 lu . 6(di) = 1.957 x 10 n and C ( d 2 ) = 5.335 x 10 1U and thus the decision is 
d = d] = (1.4.9.11). I 11 a practical implementation of (52) only the first ( .V — 1 ) symbols 
in the decoded sequence will be retained as valid symbols since the last one is actually 
the concatenation of two partially complete symbols and is likely to be in error. For our 
example above only the subsequence (1.4,9) will be retained. By properly allowing overlap 
between consecutive observation intervals, all symbols are effectively decoded. 

Approximations to (51) can be derived that further reduce complexity. An obvious 
one is to maximize with respect to d 6 A* not the sum over all m but only the largest 
term in the sum. i.e. 

v 

max V A’j., dKw .| d |. (53) 

1 — 1 

In (53) m M ( d) is the value of m — 0. 1, • • • , (Q — 1 ) that maximizes the inside sum in (51) 
• for a given d. What we have in effect is a joint estimation of the modulation sequence and 
synchronization, i.e. (53) is equivalent to 


max 


^ ^ ^ Ji ( d ) -r m • 
i=l 


(54J 


We next turn our attention to the performance achievable by the above receivers. 
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b) A bound on error probability 

As with the synchronization problem, here we are interested in finding a bound on 
the ultimate performance that can be achieved by the receivers described above. 

Our starting point is the observation that when a receiver detects a sequence at the 
wrong location, then the sequence is detected erroneously. This implies that the probability 
of sequence detection P,d is bounded by the probability of correct synchronization P ca , i.e.. 
P,d < P c , , which in turn yields 

Psd <Pr,i< 1 - (1 - 1/Q) (N+1) + ^(1 - 2/Q) iN+1) . (55) 

We observe that the PPM format imposes a severe error-floor on the sequence detection 
probability when symbol synchronization is absent. 

We now turn our attention to the symbol detection probability P sy d and prove the 
following proposition. 

Proposition 7: The symbol detection probability is bounded by 

Psyd < 1 - 0 - 1 /Q) (N+1) + ^(1 - 2 /Q) iN+1) . ( 56 ) 

Proof: Letting P ayd be the symbol detection probability in the absence of any channel 
noise, we have P syd < Psyd- We further observe that, in the absence of noise, when a 
sequence error is made (due to wrong synchronization ). all the symbols in a sequence are 
received incorrectly; otherwise, all symbols are received correctly. From this observation 
we infer that P ayd equals the probability of sequence detection in the absence of noise, 
which in turn equals P ra i.i.e.,P ay d = P r $i ■ We thus have P ayd < P r3 l which with equation 
(42) imply (56). This completes the proof. 

The symbol error probability floor predicted bv (56) in the absence of symbol syn- 
chronization can be very severe, especially for large values of Q. As an example. Q=256 is 
an alphabet size considered for some applications: if an error-probability floor of at most 
lO'" is required, a simple calculation shows that at least 435,200 slots need to be processed 
(N = 1700) in order to satisfy the performance requirements. 


V. CONCLUSIONS 


We have considered the problem of synchronization and detection of random PPM 
symbols in the presence of only slot synchronization. In characterizing the synchroniza- 
tion properties of PPM symbols imbedded in long random PPM sequences, we introduced 
the minimum distance dmtn • From our analysis we concluded that synchronization per- 
formance improves with the length of the observed sequence and degrades with increasing 
PPM alphabet size. 

In section III we derived optimal and suboptimal symbol synchronizers and a lower 
bound to the synchronization error probability; the error floor was seen to be due to the 
PPM modulation format . A way to remove the error-floor was suggested which consists of 
inserting periodically in the random data stream a pair of '‘special” symbols. The insertion 
of these symbols may be necessary to remove the error-floor, especially for large values of 
Q where it is most severe. Maximum- likelihood receivers that make sequence decisions in 
the presence of slot synchronization only are derived in section IV as well as bounds to the 
sequence and symbol detection probabilities. We observed here that the symbol detection 
probability is bounded bv the probability of correct synchronization. 

Our conclusion is that for channels where the mechanisms that can cause random slot 
shifts are such that the phase in a given observation interval cannot be reliably predicted 
from previous observation intervals, the symbol error-floor is severe, especially for large Q. 
In these cases, some signal design to eliminate the error-floor is necessary. 

VI. REFERENCES 

1. J.R. Pierce, “Optical channels: Practical Limits with Photon Counting,” IEEE Trans- 
actions on Communication, Vol. Com. -2b, No. 12, pp. 1819-1821, December 1978. 

2. D.L. Snyder and I.B. Rhodes, “Some Implications of the Cutoff-Rate Criterion for 
Coded Direct-Detection Communication Systems.” IEEE Transactions on Informat ion 
Theory. Vol. IT-26, pp. 327-338, May 1980. 

3. D.L. Snyder and C'.N. Georghiades, “Design of Coding and Modulation for a Power 


Efficient Use of a Bandlimited Optical Channel,” IEEE Transactions on Communica- 
tions, Vel. COM-31, No. 4, pp. 560-565, April 1983. 

4. C.C. Chen and C.S. Gardner. “Performance of PLL Synchronized Optical PPM Com- 
munication Systems.” IEEE Transactions on Communications, Vol. COM-34, No. 10, 
October 1986. 

5. G. Ling and R.M. Gagliardi, “Slot Synchronization in Optical PPM Communications,” 
IEEE Transactions on Communications, Vol. COM-34, No. 12, December 1986. 

6. C.N. Georghiades, “Maximum- Likelihood Symbol Synchronization for the Direct- 
Detection Optical Onoff-Keying Channel,” to appear in the IEEE Transactions on 
Communications May 1987. 

7. C.N, Georghiades, “Optimum Joint Slot and Symbol Synchronization for the Optical 
PPM Channel.” to appear in the IEEE Transactions on Communications May 1987. 

8. R. Gallager. Information Theory and Reliable Communication . Wiley tc Sons 1968. 

9. J. Lesh. J. Katz, H. Tan and D. Zwillinger, “2.5 Bit/Detected Photon Demonstration 
Program: Description. Analysis and Phase I Results,” Internal Report, Jet Propulsion 
Laboratory. 

10. J.L. Massey, “ Optimum Frame Synchronization. ’’IEEE Transactions on Communica- 
tions, Vol. COM-20, No. 2, April 1972. 

11. C.N. Georghiades and D.L. Snyder, “Locating Data Frames in Direct-Detection Op- 
tical Communications," IEEE Transactions on Communications, Vol. COM-32, No. 
2. February 1984. 

12. C.N. Georghiades, “ Joint Band and Frame Synchronization in Optical Communica- 
tions,” IEEE Transactions on Communications, Vol. COM -33. No. 4, April 1985. 


23 



APPENDIX A 


Derivation of Equation ( 1 ) 

For a given alphabet size and sequence length (N + 1) let Pr [ K\j/S\ = i] be the 
probability that K of the N symbols within the PPM binary subsequence at j are PPM 
symbols, given that the first symbol in the sequence is 5 1 = i, i = 1.2. • • • , Q.. For a given 
j, we distinguish two cases: either i < j or i > j. For i < j, the probability that K of the 
N symbols are PPM symbols is the probability that K of these symbols have pulses in the 
first j slots. Since symbols within a sequence are randomly chosen, the probability that K 
of the N symbols have pulses before or at the j — th slot is (j/Q) h ( 1 - j/Q) N ~ h ■ Since 
there are ( ^ ) sets of K symbols we have 

Pr[K-,j/S, = i] = (£)o 7Q> K <1 * = 1,2,.., j. (1.1) 

Now. for Q > i > J • in order for I\ symbol matches to occur it must be that K symbols 
have pulsed slots after the j - th slot, which implies 

PrlK-.j/S, = .] = (*)(1 -j/Q) K (j/Q) N - K . i = 0 + 11.2 ■. .<?. (2.4) 

All we need to do now is expectate over all Q first symbols i, all equiprobable. Per- 
forming the expectation yields equation (1). 
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APPENDIX B 


Derivation of Equations (36) -(41) 

We first derive equation (36). Here we are interested in the number D 0 of sequences 
for which no PPM binary subsequence is a PPM sequence. Our approach is to start 
with the Q(X — 1) possible sequences and then subtract all unwanted sequences. Using 
proposition 2, we must subtract all sequences that have all their pulses after the j — th 
slot for j = 1,2, •••,((? — 1) and all sequences that have all their pulses before the j — th 
slot, j = 1, 2, • • • , (Q — 1). It is easily seen that there are (Q — i)(W+i) of the former and as 
many of the letter; (Q — 2)^ Af+ U 0 f these sequences belong to both categories. Subtracting 
2 (Q — 1 ) <a ' ! ’ 1> from Q* A "' ,) and adding (Q — 2) ( A + to avoid subtracting sequences twice 
yields (36). (37) is obtained by dividing (36) by Q ,A + 1) . 

Let us now compute the number D ] of sequences that have exactly one PPM binary 
subsequence which is a PPM sequence at some shift j = 1. 2. • • • ,(Q — 1 ). If we let X(j) 
be the number of sequences whose PPM binary subsequence at j is a PPM sequence, we 
have 

Q - 1 

(IB) 

j=i 

For a given j, X(j) is the number of all the sequences that have either all pulses before j or 
all after j. It is easily seen that there are j N+i of the former and ( Q — j ) N+i 0 f the latter. 
From the j^' rl sequences we must subtract sequences that have all pulses after slot 1 or 
before slot (j - 1 ) since these are sequences that result in PPM sequence at some slot shift 
i < j\ there are 2 (j - l) iN + 1 - ( j - 2) A + 1 such sequences for Q > 2. Similarly, from the 
(Q - _/ ) A sequences we must subtract all those that have either all pulses after the ( j — 1 ) 
slot or all pulses before (Q— 1) ( and after j); there are [2(Q — j — l)^ -1-1 - (Q — j - 2) A +1 i 
such sequences. Denoting by h(j) the number of sequences that have exactly one match 
at j and have all pulses before j and by I>(j ) those that have exactly one match at j and 
all pulses after j . we have for Q > 2. 


h(j) 


/0, j = 1 

1 ' ~ 2 (j - 1 ) iV+1 + (j - 2) s+ \j = 2. 3, •••.«?- 1 ) 


(2 B) 


and 


r m _ JM <?--/) = ! 

2(3 ’ ~ { ( Q-j) N+1 - 2 (Q-j - 1) N+1 + (Q — j — 2 ) N+1 ,(Q - j ) = 2,3 - 1). 

(ZB) 

It is clear that X(j) = Ii(j) + h (j) and that X(j) = X(Q — j). Then 

Q - 1 Q - 1 

Bi = £xo) = 2s;/ 1 o) 

>=i >=i 

= 2 [(<? - 1)* +1 - (Q - 2)* +1 - 1] , Q > 2. 

For the special case Q = 2 it is easily seen by inspection that Dj = 2 which along 
with (4B) yields (38). Equation (39) is obtained by dividing by Q^‘ v+1) . 

Finally, to obtain (40) and (41). we observe that in order for all PPM binary subse- 
quences at j = 1.2. •••,((? — 1 ) to be valid PPM sequences it must be that all (N 4- 1) 
symbols in a sequence are the same. Since there are Q such sequences, we obtain (40) from 
which 1 41 ) follows. 


FIGURE CAPTIONS 


Figure 1: Simulation results for the synchronization probability P Ci for Q = 2, N = 
5,10- 

Figure 2: Simulation results for the synchronization probability P c> for Q — 8, N = 
20.40. 

Figure 3: Simulation results for the synchronization probability P ca for Q = 16, N = 
20.40. 

TABLE CAPTIONS 

Table 1: Comparison between simulations and the bound in (42). 
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TABLE I. COMPARISON BETWEEN SIMULATIONS 
AND THE BOUND IN (42) 




Q=2 


Q^A 


N 

SIMULATIONS 

BOUND 

SIMULATIONS 


BOUND 

0 

0.970 

0.969 

0.765 


0.773 

10 

0.999 

0.999 

0.944 


0.944 

15 

1.000 

1.000 

0.986 


0.987 

20 

1.000 

1.000 

0.996 


0.997 

25 

1.000 

1.000 

0.999 


0.999 

30 

1.000 

1.000 

1.000 


1.000 

35 

1.000 

1.000 

1.000 


1.000 

40 

1.000 

1.000 

1.000 


1.000 





0=16 


N 

SIMULATIONS 

BOUND 

SIMULATIONS 


BOUND 

O 

0.490 

0.566 

0.270 


0.447 

10 

0.741 

0.756 

0.484 


0.563 

15 

0.872 

0.870 

0.620 


0.665 

20 

0.933 

0.932 

0.720 


0.748 

25 

0.966 

0.965 

0.809 


0.813 

30 

0.982 

0.982 

0.858 


0.862 

35 

0.991 

0.991 

0.897 


0.899 

40 

0.994 

0.995 

0.925 


0.926 
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